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ABSTRACT 

Solar prominences are clouds of cool plasma levitating above the solar surface and insulated from the million- 
degree corona by magnetic fields. They form in regions of complex magnetic topology, characterized by non¬ 
potential fields, which can evolve abruptly, disintegrating the prominence and ejecting magnetized material 
into the heliosphere. However, their physics is not yet fully understood because mapping such complex mag¬ 
netic configurations and their evolution is extremely challenging, and must often be guessed by proxy from 
photometric observations. Using state-of-the-art spectro-polarimetric data, we reconstruct the structure of the 
magnetic field in a prominence. We find that prominence feet harbor helical magnetic fields connecting the 
prominence to the solar surface below. 

Subject headings: Sun: magnetic topology — Sun: chromosphere — Sun: corona — Polarization 


1. THE EMPIRICAL STUDY OE MAGNETIC EIELDS IN SOLAR 
PROMINENCES 

Solar prominences are seen as bright translucent clouds at 
~10 Mm over the solar limb because they mainly scatter light 
from the underlying disk. When seen on the solar disk, they 
appear as dark, long filamentary structures, hence called fil¬ 
aments. The main body of the filament (spine) often has 
shorter side-wards extensions (filament barbs) which, when 
seen at the limb, give the impression of extending from the 
spine to the photosphere below (prominence feet) (|Mackay 
let al.|2010| ). 

It has long been clear that a magnetic field supports the 
dense material of prominences against gravity and prevents 
them from dissipating into the faint, extremely hot corona. 
Local dips on magnetic field lines can support the plasma, 
and co uld be induced by the dense, hea vy prominence plasma 
itself ( Kippenhahn & Schliiter 1957|) , or they can exist in 


force-free ( Aulanier & Demoulin|1998[|Antiochos et al.|1994|) 
or sto chastic magnetic fields (van Ballegooijen & Cranmer 
|201Q| ). Yet, all these theoretical claims must be constrained 
by the empirical determination of magnetic fields in promi¬ 
nences. 

In the Sun, and in any astrophysical plasma in general, 
we are not able to directly measure these fields, but we are 
obliged to infer them from the light they emit. Spectro- 
polarimetry, the measurement of the polarized spectrum of 
light allows us to recover quantitative information on the 
magnetic field vector. The polarization state of observed light 
is compatible with the intrinsic (broken) symmetries of the 
emitting plasma, in particular, with the presence of a mag¬ 
netic field. Thus, for example, the emission by an isotropic 
(and therefore unmagnetized) medium is unpolarized. Polar¬ 
ized emission along a magnetic field is circularly polarized 
while, normally to the field, it is linearly polarized, as in the 
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longitudinal or transversal Zeeman effe cts, respectively (e.g. 
Land! DeglTnnocenti & Landolfi|2004| ). If light is scattered, 
additional symmetries Se broken and the dependencies of po¬ 
larization are more involved. In solar prominences and fila¬ 
ments, spectral lines are polarized by scattering a nd the Zee- 
man effect, and futher modified by the Hanle effect (|Tandberg-| 


Hans senjl 995 [ Landi^egrinnocenti & Landolfi|2004| , pro 

viding direct information on the magnetic field vector. 

Many studies have observed prominences with the aim to 
determine magnetic fields. Some maps of the magnetic field 
vector in quiescent pro minences show horizontal magnetic 
fields of ~ 10 - 20 G ( |Casini et al.||2005t |Qrozco Suarez| 
et al.|2014 ), in agreement with results obtained in the 1970’s 
and 1980^ from observations with limit ed spatial resolution 
( jSahal-Brechot et aO 19771 |Leroy|[T989] ). In contrast, verti- 
cal fields h ave been also diagnosed m prominences ( |Merenda 
et al.|2006| ). Considering prominence feet, the observation of 
vertical velocities in these structures suggest verti cal fields di¬ 
rectly connecting the spine with the photosphere ( Zirker etak] 
1998|). Also, from observations of photosphenc magnetic 
fields (not the prominence itself), the barbs are interpreted as 
a series of local horizontal dips su staining plasma at different 
heights ( |L6pez Ariste et al.|2006| ). 

In the light of these results it is clear that, from the observa¬ 
tional point of view, the precise topology of magnetic fields in 
prominences is still a matter of debate. The main reasons are 
that 1) measuring the polarized spectrum of solar prominences 
is an observational challenge, and 2) the inference of the mag¬ 
netic field vector in the Hanle regime is subject to potential 
ambiguities. In this paper, we reconstruct the topology of the 
feet of a quiescent prominence from spectro-polarimetric data 
at the He I 1083.0 nm line. We study the possible solutions to 
the inverse problem and, in contrast to previous results, we 
discard some of them using a physical constraint. We also 
propose an analytical method to be used to find the multiple 
solutions in the case of Hanle diagnostics in prominences and 
write the explicit equations in the Appendix. 


2. NEAR-INLRARED SPECTRO-POLARIMETRY AND 
MULTIWAVELENGTH IMAGING 

On April 24 2011 (10:00-13:00 UT), we performed four 
consecutive spectro-polarimetric scans of a quiescent promi¬ 
nence (located at 90E 42S) focusing at the 1083.0 nm multi- 
plet using the Tenerife Infrared Polarimeter ( |Collados||1999 ) 
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Figure 1. a, Prominence seen at the core of Hq,. This snapshot was recorded half way of the spectropolarimetric scan, b, Prominence as seen at 17.1 nm. This 
snapshot was recorded 10 h after the spectro-polarimetric scan, when the pillars of the tornadoes were more prominent (no need for sharp-marsking filter) than at 
the time of the observations, c, Reconstructed intensity map at the core of the He I 1083.0 nm line, d, e, Maps of the amplitude of linear polarization Stokes Q 
and U. The reference direction for positive Stokes Q is parallel to the limb, f, Map of the amplitude of circular polarization Stokes V. Insets show detailed views 
of the circular polarization in both feet. Insets fl-f2b: Stokes V in the region after subtracting the average value (fl and f2b) 


at the German Vacuum Tower Telescope in the Observatorio 
del Teide. We integrated 30 s per scan step to reach a polari- 
metric sensitivity of 7x 10“^ times the maximum intensity. 
Each scan of the prominenece took around 30 min. This data 
set constitute a unique time series of high polarimetric sensi¬ 
tivity and unprecedented spatial resolution (~ 470 km on the 
Sun) of a prominence. We applied standard reduction pro¬ 
cedures to the raw spectra (bias and flat-field correction, and 
polarization demodulation) to obtain maps of the four Stokes 
parameters I, Q, U, and V (Fig. Ic-f). The slit was always 
kept across the solar limb (horizontal direction in the images), 
which allo wed us to correct for seeing-ind uced cross-talk and 
stray light ( Martinez Gonzalez et al.|201^ . The seeing condi¬ 
tions were excellent, the adaptive optics system often reaching 
an apparent mirror diameter of 20 cm. This made the applied 
seeing-induced corrections to be very small, i.e., close to the 
limb, where the effects of seeing are expected to be the largest, 
the average corrections applied were 5x 10“^ times the maxi¬ 
mum intensity. 

Simultaneous images were taken with a narrow-band Lyot 
filter centered at the core of the line with a cadence of 1 
s. Thes e images were treated with blind deconvolution tech¬ 
niques ( [van Noort et al.|2005] ) to resolve very fine spatial de¬ 
tails of the temporal evolution of the prominence (Fig. la and 
the online version of the Hq, movie). Further context was pro¬ 
vided by imaging in the coronal line of Fe IX at 17.1 nm (Fig. 


lb) observed with t he Atmospheric Imaging Assembly (AIA; 
Lemen et al. 2012|) instrument o nboard the Solar Dynamics 
Observatory ( Pesnell et al.|20T^ . 


In this paper, we focus on the study of the third scan be¬ 
cause of the double-helix appearence of the prominence feet. 
Figure 1 displays the multiwavelength intensity, and polari¬ 
metric imaging of the observed prominence. Both in the He 
I 1083.0 nm scan and Hq, intensity images, one of the two 
prominence feet (f2 from Fig. If) shows a clear double-helix 
structure formed by two fibrils. A more compact, twisted he¬ 
lical structure can also be guessed in foot fl. The two promi¬ 
nence feet correspond to vertical, dark (absorption) structures 
observed at 17 .1 nm that resemble those recently nam ed so¬ 
lar tornadoes dSu et al.|2012[|Wedemeyer et al.|2013| ). They 
were observable in the AIA data for almost three more days, 
and on April 26 (23:38 UT), they suddenly erupted, showing 
a clear helical shape. Interestingly, the foot f2 presents oppo¬ 
site polarities of the Stokes V parameter at both sides of one 
fibril (f2a). This means that the magnetic field has reversed 
polarities along the line of sight at both sides of this fibril. 
Moreover, when subtracting the average circular polarization 
(i.e., the mean longitudinal magnetic field) to the other fibril 
of feet f2 (panel f2b) and to the other feet (fl), we find similar 
patterns. 

3. INFERENCE OF THE MAGNETIC AND DYNAMIC STRUCTURE 
OF THE PROMINENCE FEET 
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Figure 2. a, Inferred magnetic field topology of the prominence feet. Direc¬ 
tion of the field projected onto the plane of the sky (short lines) overplotted 
to the intensity of the He I 1083.0 nm line and magnetic field strength (b). 
Short-dashed black lines in panel a trace the axes of the fibrils in the two feet. 
In the sketches below, blue arrows represent the actual inferred field vector in 
f2a fibril. Note the opposite polarities of the magnetic field in both sides of 
the local vertical. This, and the relative inclination of the projected held with 
respect to the hbrils imply a helical held (red lines). 

We analyze the spectro-polarimetric data using the numeri¬ 
cal co de HaZeL (Hanle Zeeman Light; [Asensio Ramos et ah] 
[2008| ) to recover the full magnetic field vector and the thermo¬ 
dynamical properties of the plasma. Facing an inverse prob¬ 
lem with observational data and a large number of dimensions 
is always ill-posed. A classical inversion code such as HaZeL 
retrieves one atmospheric model that fits the observed pro¬ 
files, though others may exist. In our case, the number of 
these ambiguous solutions depend on the regime of the mag¬ 
netic field and, more importantly, on the scattering geometry. 
Our approach is to find compatible solutions of the same in¬ 
verse problem in each pixel and then select the global scenario 
physically compatible with the context. This procedure step 
is essential to reconstruct the global topology of the magnetic 
field in the prominence. 

3.1. Determination of the scattering geometry 

The angle of the emitting atom in the local vertical (or the 
observed strcuture if we assume it in the same plane) with re¬ 
spect to the line-of-sight (the scattering angle 6) is a very im¬ 
portant parameter to correctly infer the magnetic field vector 
from spectro-polarimetric signals generated from scattering 
processes. In order to determine it, we used the images of the 
Atmospheric Imaging Assembly (AIA) at 17.1 nm in which 
we identified the feet of our observed prominence as two dark 
(absorption) vertical filaments (Fig. lb). We followed these 
filaments as they entered onto the disk and detected their po¬ 
sitions (as a projection onto the plane of the sky). We measure 
the projected distance of the filaments to the limb d over time 
t and fit it with the approximate expression 

inferring the value of tumb, the time at which the prominence 
feet were at the limb. The symbol R is the radius of the paral¬ 
lel at latitude I = 42°, and co = 0.55° h“^ is the solar angular 
velocity at that latitude. The distance d was obtained as the 



Figure 3. Artistic representation of the three-dimensional magnetic field in 
the observed. In a, blue and red lines represent arbitrary, recovered field lines 
in the prominence feet; gray lines follow field lines in the spine according 
to overall models for prominence structure and are only drawn as context 
(i.e. the twist of field lines does not represent reality). Background image: 
intensity in the He I 1083.0 nm line, b, c. Side views (see the online movie 
of the three-dimensional reconstruction). 

length of the horizontal line from the foot to the limb, which 
is a good approximation close to the limb. The inferred scat¬ 
tering angle (the angle between the line of sight and the local 
vertical) is approximately given by 

0~-cosl{2dlR)ohy\ ( 2 ) 

The distance {dlR)obs obtained by substituting t = fobs, i-o., 
the time of the observations, in Eq. After this procedure, 
we obtain that the spectro-polarimetric scan was taken at a 
scattering geometry of ^ = 98°, i.e. while the prominence 
was slightly behind the limb. 

3.2. Determination of multiple solutions 

We follow a 2-step inversion scheme to obtain a robust con¬ 
vergence of the code. Since the magnetic field is of second 
order to the intensity profile, we first use the intensity pro¬ 
file alone to infer the thermodynamical quantities. On a sec¬ 
ond step, we fix the thermodynamical parameters and find the 
magnetic field vector using the information of the polarization 
profiles. 

As stated before, this may not be the unique solution to the 
problem. In order to capture all possible solutions we fol¬ 
low, again, a 2-step procedure. First, we sample the space of 
parameters with approximate (though appropriate) analytical 
expressions. Finally, we use these analytical solutions as ini¬ 
tial guesses for a second HaZeL inversion. This will allow 
us to refine the solutions in the general unsaturated regime to 
overcome the approximations we have made. 

In the case of an optically thin plasma, a normal Zeeman 
triplet, and a magnetic field in the saturated Hanle effect 
regime - when the Larmor frequency is much larger than the 
inverse of the characteristic time for scattering-, the geomet¬ 
ric dependencies of polarization are simply expressed through 
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the Stokes parameters as ( Casini et al.|2005| ): 


Qoc(3cos^ 6 b - l)sin^ 0 fiCOs 2 Ofi 
Uoc(3 cos^ 6b - 1) sin^ 0^ sin 20^ 


Vex cos 05 , 


(3) 


showing a dependence on the inclination 0 ^ and azimuth O 5 
of the magnetic field with respect to the line-of-sight (LOS), 
like in the Zeeman effect, and an additional dependence on the 
geometry of the scattering event through the inclination 6b of 
the magnetic field with respect to the local vertical (LV). The 
linearly polarized components (Stokes Q and U) are domi¬ 
nated by scattering polarization and the Hanle effect, while 
the circularly polarized component (Stokes V) is generated 
by the longitudinal Zeeman effect, and defines the polarity of 
the magnetic fi eld along the LOS ( |Landi DeglTnnocenti &| 
|Landolfi|2004b . 

Equations show a dependence of polarization on Og and 
O 5 that yielcUhe well-known 180°-ambiguity for O 5 of clas¬ 
sical Zeeman diagnostics. The additional dependence on 6b, 
characteristic of scattering processes, may yiel d two addi¬ 
tional ambiguous solutions for the magnetic field ( [Casini et ah 
|2005| ). The He I 1083.0 nm line often forms close to the sat- 
uration regime described by Eqs. [^and all (up to four) am¬ 
biguous field configurations can be dete rmined analytically 
(see the Appendix) (see also Judge |20Q7] ). 

In order to have an idea of the potential multiple solutions 
we face, we have represented in Eigurej^the so-called Hanle 
diagram: the amplitude of Stokes Q and U in terms of the 
inclination (6b) and azimuth ( 0 ^) of the magnetic field in the 
LV. The diagram is computed for a point at 15" above the 
solar limb (middle-upper parts of the observed prominence 
feet). We have assumed the saturation regime (in particular 
a field strength of 30 G), and the scattering angle determined 
by our observations (98°). The zero inclination is defined as 
a field in the LV, and the zero azimuth is defined in the LOS, 
being 05 = 180° a field directed away from the observer. 

The more vertical fields (6b < 35° and 6b > 145°) have 
four possible solutions irrespective of the azimuth value: two 
of them (the ones represented by filled and empty circles) 
are not ambiguous if circular polarization is observed, since 
it provides the sign of the longitudinal component of the mag¬ 
netic field. The other two solutions have opposite polarities 
of the field in the LV. Eields with inclinations between 35° 
and 145° have eight possible solutions: a more vertical incli¬ 
nation with the four solutions stated above and another four 
solutions with more inclined fields. In the case of our ob¬ 
servations, the sensitivity of circular polarization allows us to 
constrain the azimuth range, and hence only four ambiguous 
solutions need to be determined. The polarity of the field (in 
the LV) can not be determined. However, as we will see, this 
ambiguity is unimportant for the purposes of this paper. 


3.3. The helical magnetic field 

The ambiguities apply at each pixel of our observations. 
However, we assume that there are no tangential discontinu¬ 
ities or shocks and hence the magnetic field in the prominence 
draw continuous lines. We reconstructed four global topolo¬ 
gies of the magnetic field which could be grouped into two 
broad categories: one with fields inclined by 6b ~ 30°, 150°, 
and another one with more inclined fields (6b ~ 90°). In both 
cases, the projection of the field onto the plane of the sky is at 
an angle with the axis of each fibril that form the prominence 
feet. The fibril f2a (displayed in Eig. 1) shows a LOS polar¬ 
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Figure 4. Hanle effect diagram for the He I 1083.0 nm line computed with 
the HaZeL code assuming a 98° scattering geometry and the saturation 
regime (B=30 G). The inclination (Ob) and azimuth ( 05 ) are referred to the 
LV. The xb = ^ correspo nds to the line of sight, and xb = 0 for a radial 
vector away from the Sun JAsensio Ramos et al.|2008| see). The solid line 
represent the curves at constant inclination. I he Hanle diagram for the range 
of inclinations between 90° and 180° is the same but with a reversed color 
palette. The orange dots represent the observed polarization amplitudes.All 
the observed points have ambiguous vertical and horizontal solutions. How¬ 
ever, note that the sign of the circular polarization is only compatible with 
one of the solutions represented by hlled or empty dots. 
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Figure 5. a, LOS velocity field from the inversion of the He 1 1083.0 nm ob¬ 
servations with intensity contours (white lines) overplotted for reference, b, 
Space-time variations of Hq, intensity for an artificial slit at 5.8 Mm, parallel 
to the limb across the helical structure (gray line). Brackets mark the time 
interval during the scan of f2. Dotted yellow lines trace the periodic move¬ 
ment of the two fibrils in the double-helix. The composition of both motions 
indicate rotation of the structure. 


ity reversal at opposite sides of its axis. Similar LOS polar¬ 
ity inversions appear across the axes of the other fibrils when 
the average value of Stokes V in the region, corresponding to 
the mean LOS component of the magnetic field, is subtracted 
(Eigs. Ifl and lf2b). Put together, this points to an helical 
global topology in both families of solutions, with the more 
horizontal configurations showing a more twisted field. 
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In order to disambiguate the problem, we introduce a phys¬ 
ical constraint through a sta bility analysis. Accord ing to the 
Kruskal-Shafranov criterion ( [Hood & Priest|1979] ) a kink in¬ 
stability develops when the amount of magnetic twist exceeds 
a critical value, so that a structure is stable if 


L Be - ' 


(4) 


where r and L are the radius and the length of the fibril, 
respectively, and and Bq are the vertical and azimuthal 
components of the magnetic field, respectively. Estimating 
r = 1.7 Mm and L = 11.6 Mm for our prominence, the stabil¬ 
ity criterion yields 0.09 for the solution with the fields around 
90° and 1.28 for the fields with inclinations of 30°, 150°. 

Figure 2 shows the magnetic field strength and topology of 
the magnetic field. The magnetic field strength has very sim¬ 
ilar values throughout all families of solutions. In most of 
the prominence the magnetic field strength is below 20 G, but 
shows filamentary structures, parallel to the helix structure, of 
higher field strength (~ 60 G; Fig. 2b). The projection of the 
field onto the plane of the sky is always roughly perpendic¬ 
ular to the solar limb, and at an angle (~ 30°-50°) with the 
axis of the fibrils that form the prominence feet (Fig. 2a). All 
fibrils present magnetic field polarity reversals at both sides 
of their axis. This implies a helical field along the fibrils with 
the axis in the plane of the sky in f2a (since the average lon¬ 
gitudinal magnetic field is close to zero), and slightly tilted 
relative to the LV in the other fibrils (since polarity reversals 
are only seen when the mean longitudinal magnetic field is 
subtracted). The helical magnetic fields of the two fibrils are 
interlaced to form a double-helix that constitutes one foot of 
the prominence (see Fig. 3 and the online movie of the artistic 
representation of the field topology). The feet magnetically 
connect the spine with lower layers, in contrast to previous 
works sugges ting that feet are just a col lection of dips at dif¬ 
ferent heights [L6pez Ariste et al.|pOQ6| ). 


3.4. Motions of the prominence material 

The intensity profiles of the He I line carry information on 
the FOS velocity, opacity, temperature, and density of the 
structure. We recover all these parameters along with the 
magnetic topology using the HaZeL code. They have the 
same values for all families of solutions of the magnetic field. 

Figure 5a shows the inferred LOS velocities. Interestingly, 
they have opposite signs at both sides of both feet of the 
promine nce (fl and f2), with value s up to 2 km s“^ (similar 
values as Orozco Suarez et al.|2012 ). Assuming that the neu¬ 
tral atoms of He still trace the field lines, this pattern of veloc¬ 
ities could be a natural consequence if the prominence mate¬ 
rial is fiowing along the inferred double-helix structure of the 
prominence feet. Note that f 1 has a positive-negative Doppler 
pattern at lower heights, while it has a negative-positive one 
higher in the feet, as we would expect for a more compact 
double-helix, which implies a larger twist. 

If the material is fiowing along the helical field lines of the 
fibrils of each feet, a Doppler pattern similar to that of large 
scales should be observed at smaller scales. The upper parts 
of f2a show a positive-negative pattern. The rest of the fib¬ 
rils have a continuous increase of the Doppler velocity as we 
move from left to right across the structure. This could be in¬ 
terpreted as a negative-positive pattern if the fibrils are not at 
a 90° scattering geometry, which is very likely. 

The two fibrils of f2 exhibit a periodic motion in the Hq, 
space-time diagram, with a period of ~ 60 min (Fig. 5b), 


consistent with the periods reported in |Su et ar] ( 2012| ) for so¬ 
lar tornadoes. Assuming that the change m Hq, brightness is 
only due to plasma movements, the observed periodic motion 
could be also interpreted in terms of plasma fiowing along 
helical field lines. Rotation of the magnetic structure is very 
unlikely since the observed ~ 50 min period would imply a 
full turn of the magnetic field in less than one hour. This will 
considerably increase the twist of the structure and will soon 
make it unstable under magnetohydrodynamical instabilities 
(easily within a day). The plane-of-the-sky periodic motions 
show maximum tangential velocities of ~ 9 km s“\ which 
are far larger than the ones inferred from the Doppler effect in 
the He i line. We must be careful to directly relate those two 
velocities since Doppler measurements are generated only by 
plasma motions while changes in H^, brightness can not only 
be assigned to plasma movements but to changes in the ther¬ 
mal conditions. 


4. DISCUSSION 


The feet of the observed prominence harbor helical mag¬ 
netic fields. After assuming a simple stability criterion, we 
have a preference for the more vertical solution and conclude 
that the magnetic fields in prominence feet connect the spine 
with the underlying atmosphere. These results are in contrast 
to the scenario in which these structures are formed by a se¬ 
ries of local horizontal dips that susta i n the plasma at differ - 
ent heights ( [Lopez Ariste et al.]|20Q6[ [Aulanier et al.||1999) . 
Assuming a uniformly twisted straight cylinder, the magnetic 
field displays a twist (number of turns over its length L) of 
1.32 at each fibril. We speculate that the connectivity of the 
prominence spine, with a well-defined helicity, and the pho- 
tospheric magnetic field below, with a fluctuating topology, 
may naturally yield the kind of helical structures we find. 

The He I Doppler velocities display opposite velocities 
along the LOS at both sides of prominence feet. The H^, inten¬ 
sity displays periodic motions in the plane-of-the-sky. Using 
only the Doppler velocities or the plane-of-the-sky motions 
alone do not allow to reach any conclusion on the plasma 
motions and can lead to controversies in the literature (e.g. 


Panasenco et al. 
the He i and the 


2014[ |^ozco Suarez et al.|2012 l. If we put 
Hq, information together, we could interpret 


these observations as the material of the prominence fiowing 
along helical field lines. However, the LOS Doppler veloci¬ 
ties inferred from the He i line and the tangential velocities ex¬ 
pected for the observed period of the Ha images do not match. 
It could be that neutral H and He clouds have different ther¬ 
modynamical properties or it could be that the changes in the 
Ha brightness have an important contribution from changing 
thermal conditions, but we need more data to really under¬ 
stand this issue. 

The prominence remains stable for several days, and we 
think that the magnetic topology of the prominence feet re¬ 
ported in this paper plays a fundamental role on the stability 
of the prominence as a whole. 
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ARRENDIX 

A. AMBIGUITIES IN THE HANLE EFFECT IN THE SATURATION REGIME 

In the saturation regime of the Hanle effect, Stokes Q and U are insensitive to the field strength, but are sensitive to the geometry 
of the field. Eor a two-level atom with a = 0 ^ Jiow = 1 transition, the optically thin limit, and the saturation regime for the 
Hanle effect, the linear polarization can be written as: 

Q=^(3 cos^ 6b - l) sin^ 0 ^ cos 2 O 5 

t/=|(3cos^6>B-l)sin^©Bsin2(DB. (Al) 

The inclination and azimuth of the magnetic field in the LV are represented by the symbols 0^ and Og, respectively. In the LOS, 
the inclination and azimuth are displayed as 6b and 0 ^. These equations are formally the same irrespective of the scattering angle 
6. The dependence on the scattering geometry is implicit in the amplitude in the abscence of a magnetic field q. 

The coordinates of the magnetic field vector B in the reference system of the vertical and the reference system of the LOS are: 

B = B (sin 6b cos 0^1 -f sin 6b sin 05 j -f cos 

B = B (sin 05 cos Ogi' -f sin 0^ sin Ogj' -f cos O^k'), (A2) 

where the unit vectors are related by a simple rotation: 

i' = cos 6i - sin ^k 

k' = sin 6i -f cos ^k. (A3) 

Given that the magnetic field vector must be the same in both reference systems, we find that the following relations apply: 

sin 6b cos (pB = sin 0 ^ cos O 5 cos 6 -f cos Og -f sin 6 
sin 6b sin 0 ^ = sin 0 ^ sin O 5 

cos 6b = cos 05 cos 6 - sin 05 cos O 5 sin 6. 
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Solving the previous three equations in the two directions, we find the following transformations between the angles in the vertical 
reference system and the LOS reference system: 


and 


cos = cos 6 cos 6b + sin 6 sin 6b cos c 

sin 05 = + 


cos O 5 = 


sin O 5 = 


cos 6 sin 6b cos (pB - cos 6b sin 6 


sin 05 


sin ^5 sin 05 
sin 05 


(A5) 


cos 6b = cos 6 cos 05 - sin ^ sin 05 cos O 5 

sin ^5 = + Vl “ cos^ 6b 

cos 6 sin 05 cos O 5 + cos 05 sin 6 


cos 05 = 


sin ^5 


sin (pB = 


sin 05 sin O 5 
sin 6b 


(A 6 ) 


Note that, since 05 e [ 0 , 7 r], we can safely use the square root and take the positive value. In order to transform from one 
reference system to the other, we can compute the inclination easily by inverting the sinus or the cosinus. However, the situation 
is different for the azimuth, because the range of variation is Therefore, one has to compute the cosinus and the sinus 

separately and the decide which is the correct quadrant fo the angle in terms of the signs of both quantities. 


Four multiple solutions can exist for the Stokes Q and U parameters. The idea is that O 5 can be modified and still obtain the 
same Q and U by properly adjusting the value of 05 . It is clear that, given that the term that can be used to compensate for the 
change in the azimuth on the LOS reference system is the same for Stokes Q and U, we can only compensate for changes in the 
sign. Therefore, we have the following potential ambiguities: 

0g = 05 

= 05 - 71 jl 

Og = O 5 + 7r/2 

^'^ = ^b+7i. (A7) 

We have to compute the value of 0^ that keeps the value of Q and U unchanged. Therefore, once we find a solution to 
the inversion problem in the form of the pair (^ 5 , 05 ), we can find the remaining solutions in the saturation regime following 
the recipes that we present now. Remember that, unless one knows the sign of cos 05 (given by the observation of circular 
polarization), the number of potential ambiguous solutions is 8 . If the polarity of the field is known, the number is typically 
reduced to 4 (or 2 if no 90° ambiguity is present). 


Under this change, we have that 


A.l. = 05 


cos 20g = cos 2 O 5 , (A 8 ) 

sin 20 g = sin 2 O 5 , 
cos Og = cos O 5 , 
sin Og = sin O 5 . 

Making use of the previous relations between the angles wrt to the vertical and the LOS, we have to solve the following equation: 

(3 cos^ ^5 “ 1 ) sin^ 05 = (3 cos^ ^5 - l) sin^ © 5 , (A9) 

which can be written as: 

1^3 (cos 0^ cos ^ - sin ^ sin 0^ cos 05 )^ - 1 j sin^ 0^ (AlO) 

= 1^3 (cos 05 cos ^ - sin ^ sin ©5 cos 05 )^ - 1 j sin^ © 5 . 

After some algebra and doing the substitution t = sin 0^, we end up with the following equation to be solved: 


+ Ct^ Vl -^2 = K, 


(All) 
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where 

A = -3 cos^ 6 3 sin^ 6 cos^ O 5 

B = 3 cos^ 6-1 
C = -6 cos 6 sin 6 cos Og 

[3 (cos 05 cosO - sin ^ sin 05 cos 05 )^ - ij sin^ 05 . (A12) 

The previous equation can be solved if we make the change of variables t = ± Vz, resulting in: 

(C^ + A^)Z^ + (-C^ + 2AB)Z^ + (-2AK + B^)Z^ (A13) 

-2BKZ + K^ = 0. 

This polynomial of 4-th order can have four different solutions. From these solutions, we have to take only the real solutions 
which are larger than 0 , given the range of variation of 05 : 

reR, 0<t<l. (A14) 

Once the solutions for t are found, we make 0^ = arcsin t. Note that, for a fixed value of t, two values of 0^ are possible. We 
choose the correct one by evaluating the expressions for Q and U and testing which of the two possible choices give the values 
equal (or very similar) to the original ones. 

The angles (Ob, (Pb) are obtained by doing the transformation from (0^, O 5 ) to the vertical reference system. 

A. 2 . = O5 -r TT 

Under this change, we have: 

cos 20g = cos 2 O 5 , (A15) 

sin 20 g = sin 2 O 5 , 
cos Og = -cos O5, 
sin = - sin O 5 . 

Following the same approach, we have to solve for 0^ in 

1^3 (cos 05 cos ^ -r sin ^ sin 0 ^ cos 05 )^ - 1 j 
= 1^3 (cos 05 cos ^ - sin ^ sin 05 cos 05 )^ - 1 j 
The solution are obtained as the roots of the same equations as before but now 

A = -3 cos^ 6 + 3 sin^ 6 cos^ O 5 
B = 3 cos^ 6-1 
C = 6 cos 6 sin 6 cos O 5 

K=^3 (cos 05 cos6 - sin ^ sin 05 cos 05 )^ - ij sin^ 05 . (A17) 

The angles (6b, (Pb) are obtained by doing the transformation from (0^, O 5 -r ;r) to the vertical reference system. 

A. 3 . O5 = O5 -l- 7r/2 

Under this change, we have: 

cos 2 O 5 = - cos 2 O 5 , (A18) 

sin 2 O 5 = - sin 2 O 5 , 
cos 05 =-sin O 5 , 
sin 05 = cos O5. 

Following the same approach, we have to solve for 0^ in 

1^3 (cos 05 cos ^ -r sin ^ sin 05 sin 05 )^ - 1 j sin^ 05 (A19) 

= 1^3 (cos 05 cos ^ - sin ^ sin 05 cos 05 )^ - 1 j sin^ 05 . 

The solution are obtained as the roots of the same equations as before but now 

A = -3 cos^ 6 + 3 sin^ 6 sin^ O 5 
B = 3 cos^ 6-1 
C = 6 cos 6 sin 6 sin O 5 

= - 1^3 (cos 05 cos ^ - sin ^ sin 05 cos 05 )^ - 1 j sin^ 05 . 


sin2 0^ (A16) 

sin^ 05 . 


(A20) 
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The angles {6b, (Pb) are obtained by doing the transformation from ( 0 ^, Og + 7 r/ 2 ) to the vertical reference system. 

A.4. Og = O 5 - 7112 

Under this change, we have: 

cos 20 g = - cos 20 fi, 
sin 2 O 5 = - sin 205 , 
cos O^^sin Ofi, 
sin Og = -cos Ofi. 

Following the same approach, we have to solve for 0^ in 

1^3 (cos 05 cos ^ + sin ^ sin 0 ^ sin - 1 j sin^ 0 ^ 

= 1^3 (cos 0 fi cos ^ - sin ^ sin 0 ^ cos - 1 j sin^ 0 ^. 

The solution are obtained as the roots of the same equations as before but now 

A = -3 cos^ 6 + 2 sin^ 6 sin^ Og 
B = 2 cos^ 6 

C = -6 cos 6 sin 6 sin O 5 

= - 1^3 (cos 05 cos ^ - sin ^ sin 0 ^ cos - 1 j sin^ 05 . 

The angles (6b, (Pb) are obtained by doing the transformation from (0^, O 5 - 71 j2) to the vertical reference system. 


(A21) 


(A22) 


(A23) 



